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Squarks and gluinos have been searched for at hadron colliders in events with multiple jets and 
missing transverse energy. No excess has been observed to date, and from a comparison of experi- 
mental cross section limits and theoretical cross section predictions one can deduce lower bounds on 
the squark and gluino masses. We present an improved analysis of squark and gluino mass bounds 
which is based on state-of-the-art cross section calculations including the summation of large thresh- 
old corrections. For illustration, we consider experimental data obtained by the CDF collaboration 
at the Fermilab Tevatron and discuss the impact of the improved cross section predictions on the 
squark and gluino mass limits. 

PACS numbers; 



I. INTRODUCTION 

The search for supersymmetry (SUSY) [l|, Q is one of the central tasks at the currently operating hadron colliders 
Tevatron and LHC. In most SUSY models the squarks (g) and gluinos (g), the coloured supersymmetric particles 
(sparticles), are produced copiously in hadronic collisions. We consider the minimal supersymmetric extension of the 
Standard Model (MSSM) |, i] where, as a consequence of R-parity conservation, squarks and gluinos are produced 
in pairs and where the lightest supersymmetric particle (LSP) is stable. For large regions of parameter space the LSP 
is also neutral and weakly interacting, and provides a viable dark matter candidate. The expected MSSM signature 
at colliders is thus characterized by multiple jets from cascade decays of squarks and gluinos and by large missing 
transverse energy -ET.miss from the two LSPs at the end of each decay chain. 

Searches for squarks and gluinos at the proton-antiproton collider Tevatron with a centre-of-mass energy of \/S — 
1.96 TeV have placed lower limits on squark and gluino masses in the range of 300-400 GeV @ in a minimal 
supergravity (mSUGRA) scenario. The proton-proton collider LHC currently operating at \/S — 7 TeV has already 
significantly extended the mSUGRA mass limits to values near 800 GeV 0, [1] • The mass bounds are deduced from 
the comparison of experimental cross section exclusion limits and theoretical cross section predictions. So far, the 
interpretation of the experimental search results have relied on calculations of the squark and gluino cross sections at 
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next-to-leading order (NLO) in SUSY-QCD Recently, the accuracy of these cross section predictions has been 

improved further by the resummation of soft-gluon emission to all orders in perturbation theory |13l4l7| . The inclusion 
of soft-gluon resummation leads to a significant reduction of the unphysical renormalization and factorization scale 
dependence and, in general, increases the NLO cross sections. 

In this work we employ state-of-the-art SUSY-QCD cross section predictions including next-to-leading-logarithmic 
(NLL) soft-gluon resummation [l3l - [l7j to derive improved mass bounds on squarks and gluinos from the analysis of 
recent CDF data [i]. 

The paper is structured as follows. In section we briefly summarize the status of cross section calculations 
for squark and gluino hadroproduction and discuss the improvement obtained from soft-gluon resummation. The 
experimental CDF analysis is reviewed in section Hill The interpretation of the cross section limits in terms of squark 
and gluino mass bounds is discussed in section ITVl where we also address the treatment of theoretical uncertainties. 
We conclude in section IVl 



II. SQUARK AND GLUINO CROSS SECTIONS 



In the MSSM with R-parity conservation, squarks and gluinos are produced in pairs, 

pp/pp qq,m,qg,gg + X . (i) 

In Eq. ([T]) and throughout the rest of this paper we suppress the chiralities of the squarks q = {q_L,q_R) and do not 
explicitly state the charge-conjugated processes. We include squarks q of any flavour except for top squarks. The 
production of top squarks [1^ needs to be considered separately because of potentially large mixing effects and mass 
splitting in the stop sector jl8j |. 

The cross sections for the squark and gluino pair-production processes ([T]) have been known at next-to-leading 
order in SUSY-QCD [ol-fl^ for some time. Electroweak corrections to the 0(a1) tree-level production [l9l - |23j and 
the electroweak Born production channels of 0{aas) and 0{a^) 24, 25] are significant for the pair production of 
SU(2)-doublet squarks q^ and at large invariant masses in general, but they are moderate for total cross sections 
summed over all squark species. 

The NLO SUSY-QCD corrections to squark and gluino hadroproduction reduce the renormalization- and 
factorization-scale dependence of the predictions. In general the NLO contributions also significantly increase the 
cross section with respect to the leading-order (LO) approximation [2^-[2^ if the renormalization and factorization 
scales are chosen close to the average mass of the pair-produced sparticles. A significant part of the NLO corrections 
can be attributed to the threshold region, where the partonic centre-of-mass energy is close to the kinematic produc- 
tion threshold. In this region the NLO corrections are dominated by the contributions due to soft gluon emission 
off the coloured particles in the initial and final state and by the Coulomb corrections due to the exchange of gluons 
between the massive sparticles in the final state. The soft-gluon corrections can be taken into account to all orders in 
perturbation theory by means of threshold resummation techniques. 

Recently, such threshold resummation has been performed for all squark and gluino production processes in Eq. ([1]) 
at NLL accuracy [Tsl - flGj . For squark-antisquark production, NLL resummation has also been addressed in the 
framework of effective field theories [2^, [s^l, and the dominant next-to-next-to- leading order (NNLO) correction 
coming from the resummed cross section at next-to-next-to-leading-logarithmic (NNLL) level has been derived in 
Ref. [31I [3^ . Moreover, a formalism allowing for the resummation of soft and Coulomb gluons in the production 
of coloured sparticles has been presented in Refs. [H, [s^, and bound state effects have been studied for gluino- pair 
production in Ref. [33*1 . 

In this work we will employ the state-of-the-art SUSY-QCD predictions at NLO-fNLL for squark and gluino 
production to derive mass bounds from the comparison with experimental exclusion limits. In the remainder of 
this section we briefiy review the calculation of the NLO-I-NLL cross section presented in [Tsl - flTj and present a few 
illustrative numerical results. 

The hadronic threshold for inclusive production of two final-state particles with masses and m; corresponds to 
a hadronic center-of-mass energy squared that is equal to 5 = (mfe -I- to/)^. Thus we define the threshold variable p, 
measuring the distance from threshold in terms of energy fraction, as 

p = (2) 

s 



Our results are based on the following expression for the NLL-resummed cross section matched to the exact NLO 
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calculation of Refs. |9l-lll] 

(NLO+NLL)/ r 21 2\ (NLO) / r 2i 2\ 

^ f dNp-"" f,/,,,{N + l,fi^)f,^^,{N + l,fi^) 



27ri 



X [agl, (iV, {m^}, - agl, (TV, {m^}, ^?) ] , (3) 



where the last term in the square brackets denotes the NLL resummed expression expanded to NLO. The initial state 
hadrons are denoted generically as hi and /12, while /i is the common renormalization and factorization scale. The 
resummation is performed after taking a Mellin transform (indicated by a tilde) of the cross section, 

c'hih-,^ki{N,{m^}) = I dp p^^'^ (Th-^h2^ki{p,{'^^}) ■ (4) 
Jo 

To evaluate the contour CT of the inverse Mellin transform in Eq. ([3]) we adopt the "minimal prescription" of Ref . [sj] ■ 
The NLL resummed cross section in Eq. ^ reads 

I 

where we have introduced the hard scale — {ruk + mi)^. The colour-decomposed leading-order cross sections in 
Mellin-moment space are denoted by GfjL,ki i with / labeling the different possible colour structures. The expressions 
for these, both in moment and in momentum space, can be found in [13, [3 • The perturbative functions Cij^kij 
contain information about hard contributions beyond leading order. This information is only relevant beyond NLL 
accuracy and therefore we keep Cij^kij = 1 in our calculations. The functions and Aj sum the effects of the 
(soft-)collinear radiation from the incoming partons. They are process-independent and do not depend on the colour 
structures. They contain the leading logarithmic dependence, as well as part of the subleading logarithmic behaviour, 
and are listed e.g. in Ref. [14i] . The resummation of the soft-gluon contributions, which does depend on the colour 
structures in which the final state SUSY particle pairs can be produced, contributes at the NLL level and is summarized 
by the factor 

Ail., {QKNplp') = exp [ f^"" I ^ I?z ] . (6) 

The one- loop coefhcients Di are derived and listed in HH • 

We now present a few selected numerical results to illustrate the impact of the NLL resummation on squark and 
gluino production at the Tevatron. We compare the LO, NLO and NLO+NLL matched results and discuss the 
theoretical uncertainty due to the choice of renormalization and factorization scales, and the parametrization of the 
parton distribution functions (PDFs). The NLO cross sections have been calculated in Refs. [9-12] and are available 
in the form of the public computer code Prospino[3^. The MS-scheme with five active flavours is used to define 
the QCD coupling as and the parton distribution functions at NLO. The masses of the squarks and gluinos are 
renormalized in the on-shell scheme, and the SUSY particles are decoupled from the running of ctg and the PDFs. 

We first discuss the scale dependence of the SUSY-QCD cross-section prediction for the four different production 
processes pp — >■ qq, qq, qg, gg+X at the Tevatron. For convenience we define the average mass of the final-state sparticle 
pair m — {mi; + m;)/2, which reduces to the squark and gluino mass for qq, qq and gg final states, respectively. The 
renormalization and factorization scales are taken to be equal. Figure [1] compares the scale dependence in LO, NLO 
and NLO-I-NLL for = rrig = m = 500 GeV. The scale p is varied around the standard scale choice po = m from 
p ~ Mo/10 up to p — 5pQ. Note that the LO predictions are obtained with LO PDFs and the corresponding LO 
values for ag. We observe the anticipated strong reduction of the scale dependence when going from LO to NLO, and 
a further significant improvement when the resummation of threshold logarithms is included, in particular for gg and 
qg production. 

At the standard choice of scale p — pQ — m the cross-section predictions are in general enhanced by soft-gluon 
resummation. The NLL if- factor i^NLL = ctnlo+nll/cnlo at the Tevatron is displayed in Figure [5] for squark and 
gluino masses in the range between 200 GeV and 600 GeV. We show results for equal squark and gluino masses, 
if- factors for different mass ratios mg/rriq can be found in Ref. [l5| . The soft-gluon corrections are moderate for qq 
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FIG. 1: The scale dependence of the LO, NLO and NLO+NLL cross sections for the four different squark and gluino production 
processes at the Tevatron. The squark and gluino masses have been set to — nig = m = 500 GeV. The MSTW-2008 [sH] 
PDF has been adopted. 



production, but very significant for gg , qq and qg final states, respectively. Because of the increasing importance 
of the threshold region, the corrections become in general larger for increasing sparticle masses. The large effect 
of soft-gluon resummation for qg and gg production can be mostly attributed to the importance of gluon initial 
states for these processes. Furthermore, the presence of gluinos in the final state results in an enhancement of the 
NLL contributions since the Casimir invariants which enter the resummation formulae are larger than for processes 
involving only squarks. The substantial value of i^NLL for qq production at the Tevatron is a consequence of the 
behaviour of the corresponding NLO corrections, which strongly decrease with increasing squark mass 

Representative values for the LO, NLO and NLO+NLL cross sections are collected in Table HI The CDF analysis 
was based on cross section predictions obtained with the CTEQ6.1 J8J PDF set. In Table |T] we thus compare 
predictions for the old CTEQ6.1 PDF and for the more recent CTEQ6.6^ PDF, which we use to derive the 
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FIG. 2; The NLL Jf-factor ATnll ~ f^NLO+NLL/cNLO for squark and gluino pair-production processes at the Tevatron as a 
function of the average sparticle mass m. The MSTW-2008 [3y] PDF has been adopted. 
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TABLE I: The LO, NLO and NLO+NLL cross sections for inclusive sparticle pair production pp — )■ qq + qq + qg ^ gg + X aX 
the Tevatron (VS=1.96 TeV), including errors due to scale variation (Aa^) in the range m/2 < fJ, < 2m. Results are shown 
for two PDF parametrizations (CTEQ6.1 [11], CTEQ6.6 [11]) with the corresponding 68% C.L. PDF error estimates. 

improved cross section limits presented in this paper. 

As discussed before, we observe an increase of the cross-section prediction near the central scale when going from 
LO to NLO and a further enhancement when NLL threshold resummation is included. The scale dependence in the 
range m/2 < /x < 2m is reduced from about ±25% at NLO to about ±15% at NLO+NLL. The estimated PDF 
uncertainty is approximately 15%. The cross section is reduced when going from the old CTEQ6.1 [s^ to the more 
modern CTEQ6.6 [39^ PDF set. 

III. EXPERIMENTAL SEARCHES FOR SQUARKS AND GLUINOS 

Experimental searches for inclusive squark and gluino production are based on the study of events with multiple 
jets of hadrons and large -ET.miss in the final state. In the CDF study |6'| at the Tevatron, the SUSY samples were 
generated using the ISASUGRA implementation in PYTHIA [i^ and normalized according to NLO cross sections as 
determined using Prospino 2.0 with CTEQ6.1M PDFs and the renormalization and factorization scales set equal to 
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the average mass of the sparticles produced in the hard interaction. The uncertainties in the NLO predictions due to 
PDFs varied between 10% and 25% across the region of the squark/gluino mass plane considered, and variations by 
a factor two of the renormahzation and factorization scales changed the theoretical predictions by 20% to 25%. Final 
states with different inclusive jet multiplicity were considered to optimize the sensitivity across the quark/gluino mass 
plane. In the scenario with squark masses significantly larger than the giuino mass at least four jets in the final state 
are expected, while for a giuino mass much larger than the squark masses dijet configurations dominate. In the CDF 
analysis, three parallel analyses in event topologies with large i?T,iniss and at least four, three, and two jets in the final 
state were carried out. In each case the event selection criteria were optimized such as to maximize the sensitivity to 
the SUSY signal. As an example, for degenerate squarks and gluinos with masses around 400 GeV, S/Vb = 6 ^ was 
obtained, with a signal efficiency, defined as the fraction of signal events passing the selection criteria, of about 12%. 
Additional details on the CDF data analysis can be found in [4l|,[43l. In this paper, CDF results are revisited with 
improved NLO+NLL theoretical predictions, leading to reduced scale uncertainties, and updated CTEQ6.6 PDF sets. 

A variation in the theoretical predictions can potentially translate into a change in the signal selection efficiency if the 
mixture among subprocesses contributing to the same final state is modified significantly. The effect of updated PDFs 
and improved NLO+NLL predictions is considered separately, and the variation of the signal selection efficiency is 
studied as a function of the squark and giuino masses in three different configurations already considered by CDF [42| : 



degenerate squarks and gluinos, fixed squark mass rriq = 460 GeV, and fixed giuino mass rrig — 230 GeV. Figure [3] 
shows the results separately for each jet multiplicity. In all cases, the changes in the signal selection efficiencies are 
small, indicating that the updated theoretical predictions do not introduce a significant variation. 

The new predictions for the SUSY cross sections in the different squark/gluino mass configurations considered 
above are compared to those used by CDF in Figure |31 Again, the impact of the PDF set employed and the inclusion 
of the new NLO+NLL terms are studied separately. In general, the use of the CTEQ6.6 PDFs slightly reduces the 
theoretical cross sections, whereas the NLO+NLL calculations predict an enhanced SUSY signal cross section with 
a reduced scale dependence. As a result, the new predictions differ from the previous ones by less than 1% across 
the different squark and giuino masses considered, but present smaller uncertainties. Only in the scenario of similar 
squark and giuino masses, the new SUSY cross sections are significantly larger, and the difference varies from 0.3% 
to 10% with increasing mass between 250 GeV and 460 GeV. 



The effect of the new improved theoretical cross sections on the calculation of the 95% confidence level (C.L.) upper 
limits for SUSY production cross sections and squark and giuino masses is investigated. In a counting experiment, 
the Poisson probability of observing n events with an expected background b and a signal efficiency e is 



where L is the integrated luminosity and s is the signal cross section. Following the original CDF study, a Bayesian 
approach [43| is adopted to compute the 95% C.L. upper limits on s given the observed number of events in the 
data. The calculation includes statistical and systematic uncertainties on e and b and their correlations, and a 6% 
uncertainty on L. In order to compute the mass bounds, the theoretical uncertainties on s are also included in the 
limit calculation via an extra uncertaity on e. For each mSUGRA point considered, observed and expected limits are 
computed separately for each of the three analyses, corresponding to different jet multiplicities in the final state, and 
the one with the best expected limit is adopted as the nominal result. 

Figure [5] shows, in the case of degenerate squarks and gluinos, the new expected and observed 95% C.L. upper 
limits on the production cross section compared to previous results. As anticipated from the observation of nearly 
invariant signal efficiencies after the inclusion of the NLL terms in the calculation, the obtained cross section upper 
limits are close to those published by CDF. The results on the new squark and giuino mass bounds are presented in 
Figure[6]for the whole squark/gluino mass plane. The impact from the new theoretical predictions on the squark and 
giuino mass limits is modest and only relevant for the case of similar squark and giuino masses, extending the CDF 
excluded region by 5 to 10 GeV. 




IV. SQUARK AND GLUINO EXPERIMENTAL BOUNDS 



e-(Ls^+'>)(Lse + b) 



(7) 



^ Here S and B denote the number of signal and SM background events, respectively. 
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FIG. 3: SUSY signal selection efficiencies for different inclusive jet multiplicities, as determined using PYTHIA MC samples 
normalized according to NLO cross sections with CTEQ6.1M fs^ and CTEQ6.6 [S^ PDFs, or using NLO+NLL cross sections 
with CTEQ6.6 PDFs. Different squark and giuino mass configurations are considered and the efficiencies are presented as a 
function of squark and giuino masses. 



V. CONCLUSIONS 



Precise theoretical predictions for sparticle cross sections are essential for the interpretation of current and future 
searches for supersymmetry at hadron colliders. Recently, the NLO-QCD cross section calculation has been fur- 
ther improved by the inclusion of threshold resummation at next-to-leading-logarithmic (NLL) accuracy. The NLL 
corrections reduce the renormalization- and factorization-scale dependence of the predictions. In general the NLL 
contributions also increase the cross section if the renormalization and factorization scales are chosen close to the 
average mass of the pair-produced sparticles. 

The impact of the new improved theoretical NLO-I-NLL calculations on the experimental bounds for squark and 
giuino production has been shown using CDF results as benchmark. Improved 95% C.L. upper limits on the SUSY 
production cross section and the squark and giuino masses are presented in the case of nearly degenerate squarks and 
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FIG. 4: NLO cross section prediction using CTEQ6.1M [sl] and CTEQ6.6 ^ PDFs, and NLO+NLL cross sections with 
CTEQ6.6 PDFs. Different squark and gluino mass configurations are considered and the predictions are presented as a 
function of squark and gluino masses. The error bars indicate the PDF and renormalization and factorization scales added in 
quadrature. Note that for display purposes the NLO and NLO+NLL CTEQ6.6 cross sections are shown at mass coordinates 
shifted by 4 GeV and 8 GeV, respectively, with respect to the NLO CTEQ6.1 cross section. 



gluinos. The difference with respect to the CDF published bounds are found to be small once the new NLO+NLL 
cross section predictions are combined with updated CTEQ6.6 PDF sets. 
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CDF Run II 
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FIG. 5: Expected and observed 95% C.L. upper limit on the SUSY production cross section in the case of degenerate squarks 
and giuinos. The new limits (using NLO+NLL and CTEQ6.6 PDFs theoretical predictions) are compared to previous CDF 
published results [^. 
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FIG. 6: Observed 95% C.L. exclusion region in the squark/gluino mass plane. The new results (dashed line) are compared to 
CDF published results (solid line). 
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